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ABSTRACT 


) 

x 

The  reentry  trajectory  of  a  bluff  zero-lift  reentry  body 
was  computed  using  an  HP-9845A  computer.  To  obtain  the 
external  shape  of  an  ablating  reentry  body  as  a  function  of 
altitude,  coefficients  for  polynomial  equations  were  deter¬ 
mined.  At  six  discrete  altitudes,  the  calculated  nose  shapes 
were  identical  to  published  information  for  the  nose  shapes 
of  a  reentry  body  with  silaceous  ablation  material.  Newtonian 
flow  and  an  exponential  atmosphere  were  assumed  throughout. 

The  computer  program  was  modified  to  include  ablation  due 
to  ice  crystal  erosion,  and  the  trajectory  was  recalculated 
for  a  model  ice  cloud  5,000  feet  thick  at  an  altitude  of 
30,000  feet.  A  total  of  four  trajectories  was  computed  for 
various  values  of  mass  loss  fraction  for  erosion  due  to  ice. 
The  differences  in  range  from  the  initial  trajectory  without 
ice  erosion  were  one  to  fifteen  feet. 
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I.  INTRODUCTION 


A.  BACKGROUND 

The  Naval  Surface  Weapons  Center,  White  Oak,  Maryland, 
is  partially  responsible  for  providing  improved  sensors  for 
balloon-carried  radiosondes  which  measure  atmospheric  proper¬ 
ties  at  various  altitudes.  One  item  of  ongoing  research  at 
White  Oak,  and  initially  the  object  of  this  thesis,  is  an 
improved  humidity  sensor. 

An  object  of  contention  within  the  strategic  missile 
community  is  whether  or  not  water  vapor,  rain  or  ice,  i.e., 
atmospheric  humidity  phenomena,  significantly  modify  the 
reentry  trajectory  of  hypersonic  zero-lift  reentry  bodies. 

In  order  to  determine  the  influence  on  the  reentry  trajectory, 
the  ability  to  measure  accurately  water  vapor,  water  droplets 
and  ice  crystals  must  exist. 

The  humidity  sensor  currently  being  used  on  radiosondes 
is  a  carbon  hygristor  element  which  has  poor  time  response 
at  high  altitudes  because  of  the  inherent  sensing  delays  in 
the  hygristor  element.  For  temperatures  of  -40°C,  which 
corresponds  to  altitudes  above  about  30,000  feet,  the  time 
lag  of  the  sensing  element  can  be  of  the  order  of  fifteen 
minutes  [Ref.  1].  Therefore,  significant  errors  in  humidity 
versus  altitude  maps  are  encountered. 

Table  1  is  a  reprint  from  Brousaides  [Ref.  1]  showing  the 
time  response  of  the  currently  used  carbon  hygristor.  Various 
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step  changes  in  relative  humidity  at  several  fixed  temperatures 
yield  the  tabulated  time  responses  to  achieve  63%  or  90%  of 
the  step  input. 


Table  1.  Carbon  Element  Response  Lags 
CReproduced  from  Brousaides  [Ref.  1],  p.  16) 


RH  Level  % 
(Initial-Final) 

Temp 

Change  in  RH  Indication,  sec 
63%  90% 

30-50 

0*C 

0.  1 

0.5 

51-85 

0.3 

1.4 

85-51 

0.4 

2.2 

51-30 

0.  3 

1.6 

31-51 

-20*C 

2.6 

14 

50-80 

2.4 

21 

79-51 

2.6 

48 

51-31 

3.8 

47 

31-50 

-30*C 

15 

61 

50-72 

11 

104 

72-50 

9 

169 

50-31 

33 

163 

31-50 

-40*C 

86 

363 

50-66 

82 

272 

65-50 

88 

578 

50-30 

260 

727 

Table  1  shows  that,  as  temperature  decreases,  that  is,  as 
altitude  increases,  the  time  response  of  the  element  becomes 
so  slow  as  to  make  the  readings  virtually  meaningless  for 
any  reasonable  rate  of  upward  balloon  motion.  In  the  case 
where  ice  is  to  be  detected,  the  time  lag  of  the  sensing 
element  is  particularly  critical. 

Ine  object  of  current  research  on  the  hygristor  element 
is  to  determine  a  heat  transfer  model  of  the  sensing  element 
so  that  the  delay  times  can  be  calculated  as  a  function  of 
temperature.  By  relating  delay  time  to  temperature  and 
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temperature  to  altitude,  a  more  accurate  map  of  humidity  versus 
altitude  can  be  obtained.  If  an  accurate  heat  transfer  model 
of  the  hygristor  element  can  be  obtained,  data  processing  can 
be  used  to  yield  accurate  measurements.  Alternatively,  an 
improved  sensor  can  be  designed  once  the  exact  cause  of  the 
delay  can  be  pinpointed. 

Even  if  humidity  can  be  accurately  measured,  determination 
of  water  droplet  or  ice  crystal  concentration  and  size  distri¬ 
bution  within  a  cloud  is  difficult.  In  other  words,  even  if 
humidity  can  be  sensed  accurately,  can  this  information  be 
related  directly  to  rain  or  ice?  Perhaps  a  different  type 
of  sensor  or  several  sensors  might  be  needed  to  determine 
humidity,  cloud  composition  and  concentration. 

In  addition,  in  order  to  design  a  useful  sensor,  criteria 
for  measurement  accuracy  are  required.  The  criteria  for 
sensor  accuracy  are  determined  by  the  sensitivity  of  the 
reentry  trajectory  to  changes  in  water  vapor,  rain  and  ice 
concentrations . 

Criteria  are  not  available  for  determining  the  change 
in  a  reentry  trajectory  as  a  function  of  atmospheric  water 
content.  Therefore,  it  became  necessary  to  start  at  the 
beginning  in  order  to  learn  enough  about  the  reentry  process 
to  make  estimates  of  how  much  the  reentry  trajectory  could 
differ  with  changes  in  humidity-related  atmospheric  conditions. 
This  thesis  estimates  the  change  in  range  due  to  clouds  with 
ice  crystals. 
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Consequently,  basic  information  concerning  hypersonic 
flow  about  bluff  bodies  was  studied.  Equations  of  motion 
for  reentry  bodies  were  formulated  in  relation  to  the  current 
problem.  The  many  factors  relating  to  reentry,  such  as 
ablation,  convective  and  radiative  heating  and  the  modifi¬ 
cations  to  ablation  because  of  chemical  interaction  with 
water  were  considered. 

B .  SCOPE 

Three  separate  atmospheric  water  phenomena  affect 
reentry  ablation.  The  three  phenomena  are  water  vapor,  rain 
and  ice  crystals.  Of  the  three  phenomena,  reentry  through 
ice  clouds  is  the  most  destructive.  Consequently,  this 
thesis  studies  the  effects  of  ice  clouds. 
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I I .  THEORY 


A.  REENTRY  DYNAMICS 
1 .  Trajectory 

The  trajectory  of  interest  is  the  ballistic  reentry 
of  a  zero-lift  vehicle  typical  of  reentry  bodies  used  as 
warhead".  A  ballistic  missile  trajectory  is  typified  by 
near-orbital  reentry  velocities  and  a  steep  initial  flight 
path  angle,  i.e.,  greater  than  five  degrees. 

All  trajectories  of  reentry  bodies  pass  through 
three  operational  phases.  The  three  operational  phases  are 
defined  as  follows: 

a.  Keplerian  phase:  The  segment  of  the  trajectory 
where  gas  dynamic  forces  are  negligible  (p  »  0) . 

b.  Intermediate  phase:  The  segment  of  the 
trajectory  where  gas  dynamic  forces  are  of  comparable 
magnitude  to  other  terms  (gravity)  in  the  equations  of 
motion. 

c.  Gas -dynamic  phase:  The  segment  of  the  trajectory 
where  aerodynamic  forces  dominate  the  equations  of  motion 
for  reentry. 

Of  the  three  phases  described  above,  the  gas  dynamic 
phase  is  the  one  in  which  environmental  factors,  such  as 
clouds  or  humidity,  influence  the  trajectory.  Therefore, 
the  gas  dynamic  phase  is  the  only  portion  of  the  reentry 
which  will  be  considered. 


The  analysis  of  the  flight  of  the  reentry  body  is  a 
particular  application  of  the  general  theory  of  the  dynamics 
of  rigid  bodies  in  three  dimensions.  The  motion  of  the  center 
of  mass  can  be  separated  from  the  motion  of  the  body  about  the 
center  of  mass.  In  this  discussion,  since  a  zero-lift  body 
is  hypothesized,  the  assumption  is  made  that  the  center-line 
of  the  body  always  is  parallel  to  the  velocity  vector. 

In  addition,  the  assumption  will  be  made  that  the 
reentry  trajectory  is  planar  rather  than  three-dimensional. 
Motion  in  the  third  dimension  is,  of  course,  possible  and 
indeed  certain  for  actual  reentry  bodies  due  to  wind, 
precession,  Coriolis  force  and  so  forth.  However,  if  the 
"particle"  can  be  described  as  one  capable  of  generating  drag 
only,  then  motion  in  only  two  dimensions  is  assured.  Conse¬ 
quently,  the  reentry  body  appears  as  a  point  particle  of 
varying  mass  subject  to  the  various  forces  applied  by  the 
reentry  environment. 

Atmospheric  properties,  particularly  the  densities 
existing  at  higher  altitudes,  are  subject  to  considerable 
uncertainty.  In  order  to  provide  a  model  of  the  density 
profile  of  the  atmosphere,  two  assumptions  are  made: 

1.  The  earth  and  its  atmosphere  are  spherically 
symmetrical . 

2.  The  atmospheric  density  varies  exponentially 
with  altitude. 

The  first  assumption  is  reasonable  for  this  analysis 
since  only  a  small  portion  of  the  atmosphere  is  traversed 
for  a  ballistic  entry.  The  second  assumption  is  based 
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on  the  simple  kinetic  theory  of  an  isothermal  gas  in  a 
uniform  gravitational  field.  These  assumptions  lead  to  the 
well-known  exponential  model  for  the  atmosphere: 

P  -  P0e'By  (1) 

where:  p  =  free  stream  atmospheric  density  (slug/ft^) 
pQ  *  assumed  mean  sea  level  atmospheric  density 
B  =•  atmospheric  density  decay  parameter  (ft’1) 

y  *  altitude  of  vehicle  above  the  surface  of 

the  earth  (ft) 

B  is  equal  to  the  mean  molecular  weight  of  the  earth's 
atmosphere  times  the  local  acceleration  due  to  gravity 
divided  by  the  product  of  the  mean  atmospheric  temperature 
and  the  universal  gas  constant.  Figure  1  is  a  reprint  from 
Loh  [Ref.  2]  and  compares  the  exponential  approximation 
with  the  ARDC  (1956)  model  of  the  atmosphere. 

2 .  Equations  of  Motion 

The  geometry  of  an  IRBM  or  ICBM  trajectory  is  shown 
in  Figure  2  [Ref.  3].  Assuming  no  friction  occurs  during 
the  flight,  the  conservation  of  energy  gives: 

1/2  mV^  -  mgr  *  1/2  mV  ^  -  mgR  ■  constant  (2) 

Equation  (2)  provides  a  unique  relationship  between 
velocity  and  altitude  provided  some  initial  condition 
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Altitude  it  >10' 


Figure  1.  Comparison  of  Exponential  Approximation 

with  ARDC  Model  of  Earth' «  Atmosphere  (1956) 
(Reproduced  from  Loh  [Ref.  2],  p.  13) 


Figure  2.  Diagram  of  ICBM  Flight  Path 

(Reproduced  from  Fuhs  [Ref.  3]) 
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or  boundary  condition  can  be  used  to  find  the  constant  term. 
Equation  (2)  also  can  be  written,  if  r  3  R,  as: 


I  tC^O2  ♦  Cr  ££)2]  -  g  *  constant  (3) 

where  the  constant  in  both  equations  (2)  and  (3)  is  the  total 
reentry  energy,  i.e.,  kinetic  plus  potential  energy,  of  the 
body. 

The  conservation  of  angular  momentum  of  a  constant 
mass  also  applies  and  can  be  written  as: 


mr2u)  =  mr2  »  mRV0(O)  =  mRV^cosA 


(4) 


Boundary  conditions  for  equations  (3)  and  (4)  are  0(0)  *  60 
and  r(0)  *  R.  Time  can  be  eliminated  from  the  equations 
using  the  relation  from  above. 


dr  .  dr  d9  .  dr  RVosA 

air  aF  at  if  ^2 


Substitution  of  equations  (4)  and  (5)  into  equation 
(3)  gives: 


dr  RVosA 
3T  —p - 


2gR 


2g  R 


rRV ,cosA> 

_ 1__ _ 


2 


(6) 


21 


Let 


2  * 


[ 


RV.cosA 

<P 


r 


gR  -i 
V  cosAJ 


(7) 


Note  that 


-RV,  cosA 

dZ  =*  - ^ -  dr 


(8) 


Substituting  equations  (7)  and  (8)  into  equation  (6)  yields 


dZ  _ 

ar 


2gR  + 


(9) 


Integrate  equation  (9)  to  obtain  9 : 


(10) 


If 


<J> 


2 


2gR  + 


UD 


then  equation  (10)  becomes 
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r 


To  evaluate  the  constant  C,  symmetry  and  the  boundary 
conditions  r  »  R,  0  »  ±0Q  and  Z  *  ZQ  can  be  used.  Thus: 


Substitution  of  the  boundary  conditions  into  equation  (12) 
yields : 

z0 

cos  (9q  +  c)  *  COS  (80  -  c)  -  y  C14) 

This  implies  that  C  =  it  and  therefore: 

cos  (0  +  it)  *  -  co s 8  =  i  (15) 

At  this  point,  refer  to  equations  (2)  and  (3),  and 
note  that  for  some  value  of  initial  velocity  the  total  energy, 
i.e.,  the  constant  term,  will  be  zero.  This  value  of  initial 
velocity  corresponds  to  the  initial  velocity  at  which  the 
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particle  kinetic  energy  is  equal  to  potential  energy.  Since 
the  constant  term  is  zero,  as  the  particle  moves,  the  constant 
term  will  remain  zero,  and  kinetic  energy  will  just  balance 
the  potential  energy  for  all  values  of  r,  even  as  r  approaches 
infinity.  Thus,  the  value  of  initial  velocity  is  the  escape 
velocity.  Hence 


2  =  v  ^ 

( p  e 


2gR 


(16) 


Substitute  equations  (7)  and  (11)  into  equation  (15)  to 

get 


cosQ 


,  gR  _ 
V,  cos A 
<P 


RV.cosA 

<t> 


r 


(17) 


Define  and  substitute  the  relations  V  =  V./V  and  r ’  *  r/R 

r  “ 

and  equation  (16)  into  equation  (17)  to  arrive  finally  at 


cos9 


2  -  (2V2cos2A)/r’ 

7l  +  4VZcosZA(VZ  -  1))1/2 


(13) 


Equation  (18)  is  used  to  calculate  a  family  of  curves 


showing  range  versus  launch  angle.  In  Figure  3,  range  is 


found  as  a  function  of  launch  angle  be  setting  r* 


1.0  in 


equation  (18)  to  find  0Q  for  various  values  of  V  between 
0.3  and  0.7  and  realizing  that 

Range  *  2R0Q 

,  7 

where  R  =  2.09  x  10  ft  for  the  radius  of  the  earth. 


Figure  3.  Range  of  a  Ballistic  Missile  as  a  Function  of 
Launch  Angle  for  Various  Values  of  V  =  V  /V 
(Reproduced  from  Fuhs  [Ref.  3])  ^  e 


A  value  of  V,  can  be  deduced  from  Figure  3  by  assuming 


that  the  missile  will  be  launched  with  the  proper  velocity 
and  launch  angle  to  achieve  the  optimum  range,  i.e.,  the  peaks 
of  the  curves. 

For  example.  Table  2  shows  the  various  characteristics 
of  typical  ballistic  missiles  in  the  inventory  today.  The  C-4 
Trident  has  a  range  of  between  4000  and  4800  nm.  Referring  to 
Figure  4,  an  optimum  range  of  4100  nautical  miles  requires  a 
launch  angle  of  about  27°  and  V  of  about  22006  ft/sec. 


Table  2.  Launch  Vehicles  and  Parameters 
(Reproduced  from  Fuhs  [Ref.  3]) 


Missile 

31.1 

Icn'Jlu 

Weight 

Payload 

P.tOQfi 

nm 

Warhead 

Yield 

IOC 

Retired 

CEP 

feet 

Stage* 

A-l  Polaris 

30. 000 

1200 

No*  60 

Oct  65 

2 

A.?  Polar) j 

S4-3I * 

1500 

2 

A-3  Polar)* 

54* 31' 

2500 

2 

C-J  Poseidon 

72*34' 

T/1RV 

2300 

1800 

C-4  Trident 

72 '34' 

70 ,000 

4000 

(4U00! 

100  XT 

19797 

?  ’..ftf 

3 

0-4  Trident 

42*45. 8' 

6700 

. 

0-5  Trident 

82*45.8' 

125,000 

MRV 

7MM?A 

350XT 

1985 

(1987) 

MX 

92" 

172,000 

I4MK12A 

6000 

1983 

t  (*!«  to 

gui^tt'c ) 

(30'  '.w  ) 

Atlit 

1965 

Titan 

HK6 

MM  I 

mil 

MHIU 

3Hr.17A 

0000 

350KT 

■ 

600 

3 

Returning  to  equation  (2),  the  constant  for  the  C-4 


Trident  can  be  found  by  substituting  and  R. 


T  V  2  -  gR  =  =  K 

L  <j>  °  m 


(20) 


K  *  j  (22006)  2  -  (-32 .2)  (2.09  x  107)  =  9.15  x  108  ft2/sec2 


Since  K  now  is  known,  equation  (2)  can  be  used  to  find  the 
velocity  at  any  point  in  the  trajectory  given  r  by  rewriting 
equation  (2)  as 


V  =  [2 (K  +  gr) ] 1/2 


(21) 


If  the  portion  of  the  trajectory  near  reentry  is 
examined  (that  is,  where  r  is  not  much  different  from  R) , 
the  earth  can  be  approximated  as  being  flat  since  y/S  1 
where  y  is  altitude.  For  this  case,  the  simple  two-dimensional 
equations  of  motion  in  Cartesian  coordinates  apply,  that  is 


EFx  =  mx  (22) 

and 

IFy  =  my  (23) 
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where  x  is  distance  measured  horizontally,  and  y  is  distance 
measured  vertically. 

Since  in  the  model  the  only  force  in  the  x-direction 
is  drag,  and  the  forces  in  the  y-direction  are  drag  and 
gravity,  equations  (22)  and  (23)  can  be  written  as 


X  = 


(24) 


y  * 


(25) 


where  Dx  and  Dy.  represent  the  drag  in  the  x  and  y  directions, 
respectively.  In  finite  difference  form,  equations  (24)  and 
(25)  are 


X 


j  +  2 


2Xj+l  *  Xj  .  Dx 
2  m 


At 


(26) 


Vi 


iVi 


+  Y, 


At* 


-g 


-  JL 
m 


(27) 


The  values  of  X-^  and  Y^  and  V(o)  are  obtained  from  initial 
conditions  at  the  point  of  interest  at  the  beginning  of 
reentry  using  equations  (18)  ,  (19)  ,  and  (21)  .  The  values  of 
X2  and  Y^  are  obtained  from 
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X2  -  Xx  *  U0At  (28) 

y2  "  Yi  +  Vt  C29) 


where 


VC0)  +  VQ2  (30) 

B.  HYPERSONIC  REENTRY 

The  character  of  the  gas  dynamic  phase  of  flight  at 
near-orbital  velocities  is  far  removed  from  that  of  flight  at 
sonic  or  subsonic  velocities.  The  vehicle  has  tremendous 
energy  and  angular  momentum,  most  of  which  must  be  released 
to  the  atmosphere.  The  two  processes  available  to  transfer 
the  energy  are  manifested  as  heating  of  the  vehicle  and 
deceleration. 

No  precise  definition  may  be  given  for  the  velocity 
at  which  supersonic  flow  becomes  hypersonic  flow  because 
the  effects  vary  with  the  vehicle  geometry.  In  addition,  the 
onset  of  these  effects  occurs  gradually  rather  than  abruptly. 
Knowledge  of  Mach  number  at  which  transition  from  supersonic 
to  hypersonic  flow  occurs  is  important.  If  the  flow  is  hyper¬ 
sonic,  certain  simplifying  assumptions  can  be  made  concerning 
the  shock  shape  and  position  relative  to  the  body,  which 
relate  directly  to  drag  and  pressure  coefficient.  These 
assumptions  are  discussed  below. 
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During  hypersonic  flight,  there  will  be  regions  of  the 
flow  fields  about  a  typically  blunt-nosed  body  where  the 
flow  decelerates  and  in  which  subsonic,  transonic  and 
supersonic  flows  all  exist  simultaneously.  For  the  purposes 
of  this  discussion,  a  Mach  number  greater  than  or  equal  to 
four  will  be  considered  hypersonic. 

Any  bluntness  of  the  nose  will  manifest  itself  as  a 
detached  shock.  A  pointed  nose  will  melt  or  ablate  due  to 
the  intense  heat  generated  at  the  point  and  therefore  will 
become  blunt.  Consequently,  the  blunted  nose  generates  the 
detached  shock.  The  most  striking  effect  of  hypersonic  flight 
is  the  formation  of  a  shock  ahead  of  virtually  all  bodies. 

The  behavior  of  the  shock  for  the  flow  past  a  blunt  body  of 
revolution  as  the  Mach  number  increases  may  be  illustrated 
as  in  Figure  4,  taken  from  Cox  and  Crabtree  [Ref.  4],  for 
flow  of  a  perfect  gas  about  a  sphere. 

As  can  be  seen,  the  shock  stands  off  from  the  body  at  all 
Mach  numbers,  and  the  shock-detachment  distance  decreases 
as  the  Mach  number  increases.  The  limiting  shock  (M  ■+  lies 
close  to  the  body  and  is  very  nearly  concentric  over  most  of 
the  surface  of  the  sphere.  Just  behind  the  shock  wave  where 
the  shock  is  nearly  perpendicular,  the  flow  is  brought  almost 
to  rest  and  is  subsonic.  As  the  flow  passes  around  the  body, 
the  flow  accelerates,  passes  through  sonic  speed  and  finally 
accelerates  to  supersonic  velocities  as  is  illustrated  in 
Figure  5,  again  taken  from  Cox  and  Crabtree  [Ref.  4], 
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Figure  4.  Shock  Shapes  for  Flow  Past  a  Sphere  at  Different 
Mach  Numbers  (y  =  1.4).  (Reproduced  from  Cox 
and  Crabtree  [Ref.  4),  p.  9) 

At  the  high  temperatures  experienced  behind  the  shock  in 
hypersonic  flow,  the  assumption  of  a  perfect  gas  cannot  be 
made  as  the  energies  of  the  gas  associated  with  temperature 
are  comparable  to  the  energies  of  various  molecular  and  atomic 
processes  such  as  ionization,  dissociation  and  vibration. 
Therefore,  characteristics  of  the  gas  such  as  the  commonly 
accepted  value  of  the  ratio  of  specific  heats  for  air,  i.e., 
y  *  1.4,  must  be  viewed  with  skepticism  and  taken  into  account 
when  calculating  the  flow  field. 

However,  the  Mach  number  can  be  associated  with  the  ratio 


of  the  kinetic  energy  and  the  random  thermal  energy  by  the 
relation 


(velocity) 2/2 


V2/2 


2 

(speed  of  sound)  /2 


T7Z  =  y^T7T 


kinetic  energy 
thermal  energy 


(31) 


Figure  5.  Subsonic  and  Supersonic  Regions  for  Flow  Past  a 
Blunt  Body.  (Reproduced  from  Cox  and  Crabtree 
[Ref.  4] ,  p.  65) 

For  large  Mach  numbers,  regardless  of  the  value  of  y,  the 
kinetic  energy  is  much  larger  than  the  thermal  energy.  This 
implies  that  the  random  motion  associated  with  collisions 
between  atoms  and  molecules  is  of  relatively  little  importance 
compared  to  the  kinetic  energy  of  the  particles  and  therefore 
the  gas  may  be  thought  of  as  a  stream  of  non-interacting 
particles.  The  stream  of  non-interacting  particles  was  first 


conceived  by  Newton  for  a  "rare”  gas  composed  of  particles  with 
large  distances  between  them.  Newton  showed  that  if  the 
particles  do  not  interact  with  one  another  and  at  large 
velocities,  then  the  resistance  of  a  body  into  which  the 
particles  collide  will  vary  as  the  square  of  the  velocity. 

In  addition,  he  calculated,  by  assuming  that  there  was  no 
transfer  of  tangential  momentum  and  that  the  normal  momentum 
of  the  particles  is  destroyed  at  the  surface,  that  the  force 
on  a  plate  at  incidence  to  the  flow  varies  as  the  square  of 
the  angle  of  incidence.  More  recently  these  assumptions  have 
been  used  to  show  that  the  pressure  on  the  plate  is  given  by 

AP  -  pV2sin2ij>  (32) 


or 


Cp  *  2sin2<l> 


(33) 


Also,  for  flow  over  wedges  at  high  Mach  number,  the 
assumptions  were  used  to  show  that  the  relation  between  the 
shock  angle  (u)  and  the  body  angle  (ip)  is 


u  ,  y  *  i 

*  ~r- 


(34) 


\  Therefore  as  y  -  1 ,  the  body  and  the  shock  coincide  and 
the  density  ratio  across  the  shock  becomes  infinite.  Thus 
there  is  a  layer  of  high  density  gas  lying  near  the  surface. 


This  is  the  "Newtonian  Shock  Layer"  or  entropy  layer  since 
very  large  changes  in  entropy  also  occur. 

Lees  [Ref.  5]  adjusted  the  Newtonian  Formula  [equation  (33)] 
empirically  so  that  the  value  at  the  stagnation  point  on  the 
nose  is  exact.  Lees'  result,  which  is  called  the  "Modified 
Newtonian  Formula,"  is  given  by 


C 


P 


cp  s in 2tp 


(35) 


where  C  is  the  stagnation  point  pressure  coefficient.  The 
po 

value  of  C  is  given  by  the  Rayleigh  formula  developed  from 
po 

the  integral  Bernoulli  equations  and  the  shock  equations  and 
is 


Y  ♦  1 
2yM  -  y 


Y  +  1 
2 


l 


(36) 
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As  in  the  Newtonian  formula,  as  M  ♦  «  and  y  1 ,  the  value 
of  C  approaches  two. 


There  are  various  other  modifications  to  the  Newtonian 
Formula  including  the  Busemann  Formula,  the  Modified 
Newtonian-Busemann  Formula  and  others  but,  as  is  shown  in 
Figures  6  and  7,  a  comparison  of  the  surface  pressure 
coefficients  for  the  flows  over  a  sphere  and  a  cylinder  at 
M  *►  ®  and  y  *  1.4  confirms  that  the  Modified  Newtonian 
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Surface  Pressure  Distributions  for  Flow  Past  a 
Circular  Cylinder  (M  ■  <*>,  y  *  1.4). 

(Reproduced  from  Cox  and  Crabtree  [Ref,  4], 


Formula  is  in  best  agreement  with  the  actual  numerical 
calculations . 

1 .  Aerodynamic  Heating  and  Ablation 

During  hypersonic  reentry,  excessive  local  heating 
is  a  serious  problem.  For  the  ballistic  vehicle,  where  the 
lift-to-drag  ratio  is  zero,  very  high  heat  fluxes  are 
encountered  but  for  a  relatively  short  duration,  i.e.,  five 
to  100  seconds.  In  practice,  the  heat  and  mass  transfer 
phenomena  associated  with  this  reentry  mode  are  nonlinear 
functions  of  time  due  to  the  nonlinear  environmental  parameters, 
chemical  reactivity,  mechanical  and  acoustical  vibration, 
gas  temperature,  mode  of  energy  transport,  total  heat  load, 
peak  heating  rate,  shape  of  the  heat  pulse,  particle  impinge¬ 
ment,  material  properties,  etc.  Of  primary  concern  is  the 
means  of  accommodating  the  energy  converted  to  heat. 

The  thermal  protection  of  the  reentry  body  can  be 
accomplished  in  three  ways:  (a)  absorption,  (b)  radiation, 
and  (c)  reduction  of  aerodynamic  heating  prior  to  reaching 
the  vehicle.  The  added  requirement  to  minimize  body  weight 
in  order  to  achieve  greater  missile  range  has  often  led  to 
the  choice  of  materials  designed  to  change  phase  when  heated 
and  to  vaporize  at  the  surface  exposed  to  the  hypersonic 
environment.  The  energy  absorption  achieved  through  such 
processes  has  proved  to  be  highly  effective. 

The  energy  transfer  mechanisms  between  the  boundary 
layer  gas  and  the  ablating  surface  are  shown  in  Figure  8, 
reprinted  from  Hurwicz,  Kratsch,  and  Rogan  [Ref.  6].  The 
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Figure  8.  Heat  Flux  Balance  at  the  Surface. 
(Reproduced  from  Hurwicz,  Kratsch,  and  Rogan  [Ref.  6],  p.  9) 


energy  balance  at  the  surface  established  the  boundary  con¬ 
dition  for  solution  of  the  boundary  layer  equations.  The  net 
heat  flux  through  the  surface  to  the  interior  of  the  ablator 
(4S)  equals  the  summation  of  the  aerodynamic  heat  flux  to 
the  surface  from  the  boundary  layer  (-q  ) ,  the  energy  trans¬ 
fer  away  from  the  surface  by  normal  gas  mass  flow  (m  h  ) ,  the 
absorption  of  incident  radiation  (aw4ra(j)  ,  the  reradiation  of 
energy  from  the  surface  (+ewaTw4) ,  the  energy  transfer  by 
mechanical  erosion  (ablation)  (^rhg),  the  energy  transfer  to 

the  surface  by  condensed  material  (m  h  ) ,  and  the  energy 

s  cs 

transfer  to  the  surface  by  internal  gas  products  (m  h  ) . 

gg  gg 

Note  that  the  aerodynamic  heat  flux  (q  )  is  the  combination 


of  convection  to  the  surface  as  well  as  chemical  energy 
transfer  due  to  surface  chemical  reactions  from  mass  transfer. 
The  preceding  factors  can  be  expressed  as  an  energy  equation 


-q 


w 


-  m  h 
w  w 


+  m 


c  c 
s  s 


+  m 


-  e  aT 
w  w 


awcirad 


-  m  h 
r  s 


(37) 


Equation  (37)  is  applicable  for  both  laminar  and 
turbulent  flow  and  couples  the  energy  balance  to  the  boundary 
layer  solution  which  is  required  to  obtain  species  concentra¬ 
tions  and  the  heat  and  mass  transfer  analysis. 

Equation  (37),  although  general,  does  assume  thermo¬ 
chemical  equilibrium  at  the  surface  and  in  the  boundary 
layer.  Such  equilibrium  may  not  exist  in  reality  since  the 
mass  ablation  rate  probably  is  reaction  rate  controlled  by 
the  ablation  material,  its  char  products,  and  the  boundary 
layer  species.  In  this  case,  gaseous  ablation  products  can 
be  determined  in  principle  through  the  simultaneous  solution 
of:  (a)  all  surface  thermochemical  equilibrium  reactions 

expressed  in  partial  pressures  and  equilibrium  constants, 

(b)  equations  of  conservation  of  elemental  species  at  the 
surface,  and  (c)  Dalton's  law  of  partial  pressures. 

Even  though  conceptually  the  surface  thermochemical 
equilibrium  reactions  can  be  expressed  explicitly,  the  actual 
products  of  the  reactions  are  difficult  to  determine  because 
of  the  elevated  temperatures  and  the  number  of  interacting 


elements.  If,  for  example,  carbon  is  used  as  an  ablator, 
the  products  normally  considered  are  C,  C2>  C^,  C^,  C^,  CO, 
CO^,  CN,  0^ ,  iV2 .  As  can  be  seen  from  Figure  9,  the  concen¬ 
trations  of  the  gases  change  as  temperature  (or  velocity) 
changes.  In  addition,  if  water  vapor  is  included,  as  is  of 
interest  in  this  discussion,  another  element,  hydrogen,  must 
be  included  as  a  possible  reactant.  Thus,  entire  organic 
chains  of  molecules  may  be  formed. 


Figure  9.  Mole  Fraction  Concentration  behind  Normal  Shock 

(Reproduced  from  Cox  and  Crabtree  [Ref.  4],  p.  10) 

In  view  of  the  complexity  of  ablation  with  non¬ 
equilibrium  chemical  compositions,  information  concerning 
nose  shape  was  sought  in  lieu  of  direct  calculation.  Informa¬ 
tion  was  found  relating  nose  shape  to  specific  altitudes.  The 
nose  tip  chosen  was  composed  of  a  silaceous  material. 
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Ablation  occurred  during  reentry  through  the  atmosphere 
including  mechanical  erosion.  Further,  turbulent  flow  was 
assumed  to  exist  from  the  start  of  reentry  and  continue 
throughout  the  flight.  Figure  10  shows  the  family  of  nose 
shapes  selected  for  this  analysis.  The  material  loss  due  to 
ablation  is  excessive  compared  to  modern  carbon  nose  tips; 
however,  the  effects  of  nose  shape  change  on  the  trajectory 
were  to  be  examined.  Therefore,  the  larger  the  ablation 
rate,  the  more  the  trajectory  could  be  expected  to  change 
from  the  non-ablating  entry. 


Figure  10.  Nose  Shape  of  Reentry  Body  at  Various  Altitudes. 

(Reproduced  from  Hurwitz,  Kratsch,  and  Rogan 
[Ref.  6] ,  p.  16) . 


Recent  experiments  have  shown  that  interesting  after 


body  ablation  patterns  are  formed  during  reentry. 


Cross-hatching  or  other  similar  spiral  ablation  patterns 
occur,  and  related  roll  torques  are  induced  producing  erratic 
perturbations  in  roll  behavior.  Such  phenomena  were  not 
considered  in  this  analysis. 


The  drag  of  a  reentry  body  experiencing  Newtonian 
flow  can  be  expressed  as: 


Drag  =  q  C  sini|>  dA 

4  p 


where  q  represents  dynamic  pressure,  Cp  represents  pressure 
coefficient,  i p  is  the  local  slope  of  the  body,  and  A  is  an 
element  of  area.  Given  the  nose  shapes  in  Figure  10  and 
assuming  that  the  shape  of  the  nose  can  be  determined  for 
every  altitude,  the  drag  at  any  altitude  can  be  computed 
using  equations  (35)  and  (36).  For  an  axisymmetric  body 
where  r  *  f(x),  the  body  can  be  divided  into  infinitesimal 
segments  of  X;  and  the  drag  can  be  found  on  each  elemental 
circular  element.  Figure  11  shows  a  typical  element  in  the 
x-r  plane.  Figure  12  shows  the  same  element  as  seen  along 
the  longitudinal  axis  and  where  R  is  the  base  radius  of  the 
reentry  body.  Therefore,  the  incremental  area  is 


2irr  dr 


pw 


Hence,  equation  (38)  can  be  written  for  the  elemental  drag 
and  becomes 


d(Drag) 


r+dr 

C  2:rrdr 
P 


(40) 


The  total  drag  on  the  body  is  then  the  sum  of  the 
elemental  drags  over  the  length. 


x=L 

Drag  *  ")  '  qCp7r(r22  -  rj2) 
x=o 


(41) 


Since  the  Newtonian  flow  model  applies  over  the 
length  of  the  nose,  the  "afterbody"  can  be  included  in  the 
nose  shape  equations  if  properly  written;  hence  equations  (38) 
to  (40)  will  define  the  drag  in  all  regions  along  the  body. 

For  projectiles  which  fly  at  low  supersonic  velocities, 
the  base  drag  is  very  important  when  considering  the  total 
drag  on  the  body.  However,  at  hypersonic  velocities,  the 
pressure  increase  across  the  shock  at  the  nose  is  so  large, 
of  the  order  of  one  hundred  atmospheres,  that  the  low  pressure 
at  the  base  of  the  reentry  body  is  comparatively  insignificant. 
In  this  analysis,  the  base  pressure  drag  was  neglected. 
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C.  HUMIDITY  EFFECTS 


As  noted  previously,  the  effects  of  water  vapor,  water 
droplets  and  ice  crystals  in  the  atmosphere  lead  to  extremely 
complex  modifications  to  all  the  equations  previously  cited. 
Considering  water  vapor  only,  density  of  the  air  is  changed. 

As  a  result,  the  modified  density  influences  reentry  body 
velocity,  ablation  rate,  etc.  Also,  as  mentioned  before,  more 
complex  chemical  equations  are  required  to  determine  ablation 
products  and,  hence,  heat  loss  and  absorption.  The  total 
effect  on  drag,  for  instance,  would  be  the  sum  of  the  effects 
on  each  parameter  affecting  drag;  for  example: 


Drag 

Drag 


density 

density 


P 


q 


A  area 
area 


(42) 


If  rain,  i.e.,  suspended  water  droplets,  also  is  included 
in  the  analysis,  mechanical  erosion  must  be  considered. 
Additionally,  when  the  droplets  cross  the  shock  wave  in  front 
of  the  vehicle,  scattering  and  annihilation  will  occur;  and 
only  some  fraction  of  the  droplet  mass  actually  will  contact 
the  surface  of  the  body.  The  interaction  of  water  vapor  and 
rain  with  the  reentry  body  is  not  considered  in  this  thesis. 

The  interaction  of  ice  in  high  altitude  clouds  would  seem 
to  be  a  worse  case  of  humidity  effect  interaction  with  a  reentry 
vehicle  since  the  assumption  can  be  made  that  most  of  the 
typical  ice  particles  will  survive  after  passing  through  the 
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shock  front.  Hence,  mechanical  erosion  should  be  increased 
greatly.  Experiments  have  shown  that  the  mass  ablated  can 
be  many  times  the  impact  mass  of  the  typical  ice  particle. 

Equations  also  are  simplified  since  the  time  in  the 
boundary  layer  is  so  short,  of  the  order  of  microseconds, 
that  the  ice  will  not  undergo  a  change  in  phase  before 
colliding  with  the  body.  Therefore,  changes  in  ice  crystal 
properties  and  chemical  interactions  were  neglected. 

Since  information  was  not  available  concerning  the 
actual  change  in  nose  shape  due  to  collisions  with  ice 
particles,  the  assumption  was  made  that  ablation  shape  was 
not  changed  but  merely  was  increased  in  rate  of  mass  loss. 
That  is,  the  identical  equations  describing  the  nose  shape 
were  used  in  both  the  baseline  trajectory  model  and  the  ice 
cloud  trajectory  model.  However,  as  mass  loss  increased 
due  to  impact,  the  nose  shape  was  changed  according  to  those 
equations  until  the  proper  mass  loss  was  achieved.  This 
assumption  is  no  doubt  incorrect  in  that  it  would  be  expected 
that  the  mass  loss  would  be  much  greater  near  the  stagnation 
point  and  much  less  along  the  perimeter  of  the  body. 


III.  SIMULATION 


A.  BASELINE  TRAJECTORY 

In  order  to  evaluate  the  consequences  of  reentry  through 
ice  clouds,  a  reference  trajectory  was  needed  to  compare  such 
parameters  as  range,  impact  velocity,  changes  in  drag  coeffi¬ 
cient,  etc.  The  equations  of  motion  were  presented  previously, 
but  initial  conditions  and  boundary  conditions  are  required 
to  solve  for  a  specific  trajectory. 

1 .  Selection  of  a  Typical  Reentry  Vehicle 

The  model  selected  as  the  reentry  vehicle  conforms 
to  the  nose  shapes  shown  in  Figure  10.  In  addition,  realistic 
values  of  ballistic  coefficient  and  drag  coefficient  were 
required.  Modern  ballistic  reentry  vehicles  have  ballistic 
coefficients  (weight  divided  by  the  quantity  drag  coefficient 
times  base  area)  in  the  range  of  1500-2000  pounds  per  square 
foot  and  drag  coefficients  slightly  greater  than  0.1. 

The  weight  of  the  reentry  vehicle  was  not  known, 
but  a  specific  gravity  of  1.4  was  selected  as  typical. 

During  the  course  of  the  initial  iterations  of  the  computer 
model,  the  body  length  was  changed  to  add  enough  mass  so  that 
the  ballistic  coefficient  was  near  the  desired  value  of  2000. 
The  body  shape  was  achieved  by  extending  the  after-section  of 
the  initial  nose  shape  shown  in  Figure  10.  Figure  15  shows 
the  final  model  shape  prior  to  reentry. 


Figure  13.  Model  Reentry  Vehicle  Shape 


2 .  Initial  Conditions 

Referring  to  Figure  4,  the  center  trajectory  was 

chosen  as  a  representative  trajectory  of  a  typical  Navy 

reentry  vehicle  which  might  be  launched  by  Polaris  or  Trident. 

This  curve  corresponds  to  a  maximum  range  of  about  4100 

nautical  miles  if  launched  at  an  angle  of  27°  with  an  initial 

velocity  of  22006  feet  per  second.  Substituting  the  values 

for  launch  angle  and  initial  velocity,  equation  (20)  yields 
8  2  2 

K  =  9.15  x  10  ft  /sec  as  shown  previously. 

Since  deceleration  and  ablation  have  not  yet  begun 
at  250,000  feet,  this  altitude  was  chosen  as  the  initial 
altitude  to  begin  reentry.  For  an  altitude  of  250,000  feet, 
the  value  of  r,  the  distance  to  the  center  of  the  earth,  is 
2.115  x  107  feet;  the  value  for  r  was  substituted  into 
equation  (21)  to  yield 

V  =  [2 (K  +  gr)]1/2  =  2 [  9.15  x  108  -  32.2(2.115  x  107)]1/2 

=  21,632  ft/sec  (43) 

At  an  altitude  of  250,000  feet  atmospheric  density  is  too  low 
to  cause  drag.  Also,  at  that  altitude  r  is  not  much  different 
from  R,  the  radius  of  the  earth.  Therefore,  the  coordinate 
system  can  be  shifted  from  r,9  to  x,y.  Since  x  forces  are 
zero,  the  x  component  of  velocity  will  remain  constant  and 
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»  fe  - 


will  be  equal  to  the  launch  velocity  times  the  cosine  of  the 
launch  angle 


V  -  V  cosA  -  22006  cos  27°  -  19,60^  ft/sec  (44) 

X  (J> 

At  250,000  feet,  the  x-component  of  velocity  must  be 
the  same  value;  and  since  the  value  of  V  is  known,  the 
reentry  angle  A  can  be  determined  from 


A 


cos 


V 

X 


cos 


-1 


19607  .  ,r 

1TFTZ  ~  “5 


04' 


(45) 


Checking  these  values  with  the  curve  selected  in 
Figure  10  and  using  equation  (18) , 


cos  9, 


[1 


1  -  2  (.  6 ) 2  cos2(27°)/l 
4(  .06) 2  cos2  (27°) ( .62 


-  D] 


i/: 


=  0.32" 


Therefore,  the  initial  angular  distance  from  apogee  is 
9q  a  34.21°.  Thus  total  range  using  equation  (19)  is 

Range  =  2R9  =  2.CL Q.9„*-i°  3  (5 =  4160  nm 

0  6000  T80 

which  corresponds  very  well  to  the  expected  value  of  4100  nm. 
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The  range  the  reentry  vehicle  is  expected  to  fly  can 


be  found  from  the  reentry  location  to  impact.  To  find  the 
range  for  a  no-drag  trajectory  recall  that  the  arc  length 
traveled  is 


9'  -  0Q  -  e  (46) 

Again,  using  equation  (18)  and  substituting  values  for 
V,  A,  and  r'  at  the  entry  point 

COS9  .  - 1 

[1  ♦  4(.6)“Cos*(25.04°) ( . 6“  -  l)]  ' 

and  9  =  cos'1f.S43)  =  32.54°. 

The  angular  arc  of  flight  using  equation  (46)  is 

9’  -  9q  -  9  =  34.21°  -  32.54°  »  1.67° 


and  using  equation  (19) 


Since  the  sequence  of  equations  above  assumed  no  drag,  the 
range  should  be  the  maximum  expected  range  of  the  model 
reentry  vehicle. 


Table  3 

Summary  of  Initial  Conditions  at 


an  Altitude  of 

250,000  Feet 

Initial 

Velocity 

21,632 

f t/sec' 

Initial 

Entry  Angle 

2 

5° 

Initial 

x  velocity 

19,607 

f t/seci 

Initial 

y  velocity 

9,138 

ft/sec 

Maximum 

range 

101.5 

nm 

— > 

3 .  Nose  Shape 

The  nose  shape  is  the  primary  factor  in  determining 
total  vehicle  drag.  Since  the  nose  shape  was  given  only  for 
discrete  values  of  altitude  as  in  Figure  10,  the  nose  shape 
was  constructed  artificially  as  a  continuous  function  of 
altitude.  This  was  done  by  using  the  HP-9845A  computer  along 
with  the  HP-9874A  Digitizer.  Each  of  the  six  nose  shapes  was 
divided  into  two  segments;  one  segment  was  from  the  nose  tip 
to  2.4  inches  aft  of  the  nose  tip,  and  one  segment  was  from 
2.4  inches  aft  of  the  nose  tip  to  about  20  inches  aft  of  the 
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nose  tip.  Another  function  was  required  to  describe  the 
recession  of  the  nose  tip. 

Each  nose  shape  was  traced  with  the  digitizer 
inputting  data  points  into  the  computer.  A  curve  fitting 
routine  was  used  to  construct  curves  of  nose  radius  as  a 
function  of  length  along  the  centerline  of  the  reentry  body 
using  the  data  from  the  digitizer.  The  best  fit  curve  for 
the  segments  from  the  nose  tip  to  2.4  inches  aft  of  the  nose 
tip  was  a  fifth-order  polynomial  with  constant  coefficients. 
The  segment  from  2.4  inches  aft  of  the  nose  tip  was  described 
by  a  straight  line.  Appendix  A  contains  graphs  of  each  of 
the  six  nose  shapes  (two  segments  each)  comparing  the  data 
input  with  the  plot  of  the  equations  which  describe  the  data. 
As  can  be  seen,  the  graphs  very  closely  approximate  the 
actual  nose  shapes. 

The  nose  shapes  thus  were  defined  mathematically  as 


1 

! 


5 

r  =  2  A,  x1 

i=o  1 


(47) 


for  0  <  x  <  2.4  inches  and  as 


r 


Z  B  •  x1 
i=o 


(48) 
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for  x  >  2.4  inches  where  r  is  the  local  radius  of  the 
reentry  body  and  x  is  the  distance  along  the  centerline. 

Using  the  coefficients  calculated  in  equations  (47) 
and  (48)  ,  plots  of  each  coefficient  versus  altitude  were 
constructed.  For  example,  all  the  fifth-order  coefficients 
of  equation  (47)  were  plotted  at  the  altitudes  applicable 
to  the  particular  nose  shapes.  The  coefficients  in  equation 
(47)  thus  were  defined  by 


A. 


l 


6 

Z 

j=o 


Cijh- 


(49) 


where  h  is  altitude  in  feet.  The  coefficients  in  equation 
(48)  were  defined  by 


B. 

l 


6 

Z 

j=o 


D.  -hJ 
ij 


(50) 


Division  of  the  coefficient  curves  into  two  portions  was 
necessary;  equations  (49)  and  (SO),  sixth-order  polynomials, 
applied  below  90,000  feet.  Equations  (51)  and  (52) 


A. 

i 


1 

Z 

j=0 


E.  .h 
ij 


J 


(51) 


54 


(52) 


mm 


B. 

i 


\  FijhJ 

J  =o 


were  used  above  90,000  feet. 

Equations  (47)  through  (52)  can  be  combined  to 
yield  the  nose  shape,  the  radius,  as  a  function  of  altitude 
in  the  various  regions  of  applicability.  Table  4  shows  the 
combined  equations  (53)  through  (56)  and  the  regions  of 
applicability. 


Table  4 

Summary  of  Nose  Shape  Equations 


Altitude 


Equal  to 
|  or  Below 
]  90,000 
i  feet 


I  Above 

;  90,000  :  r 

I  feet 


0  <  x  <  2.4  (inches) 


x  >  2.4  (inches) 


5  6 


Z  Z  C-.h^x1  (53)  |  r 
i=0  j=0  ^ 


1  6 
i=0  j  =0 


Z  Z  D- -h^x1  (54) 


Z  Z  E..hV  (55) 
i=0  j  =0  ^ 


1  1 


r  = 


Z  Z  F-.hV  (56) 
i=0  j  =0 


Appendix  B  compares  graphs  of  the  polynomials 
describing  the  A^  coefficients  versus  altitude  with  corre¬ 
sponding  graphs  of  the  digitizer  data  input  used.  Appendix  C 
contains  similar  graphs  of  the  nose  tip  recession  versus 
altitude.  All  graphs  show  that  the  equations  adequately 
represent  the  described  data. 
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The  result  of  the  above  analysis  was  to  allow  the 
insertion  of  any  altitude  into  equations  (S3)  through  (56) , 
yielding  the  corresponding  nose  radius  as  a  function  of  the 
centerline  distance  from  the  nose  tip. 

4.  Flow  Chart 

Using  the  initial  entry  conditions  and  nose  shape 
versus  altitude  information  previously  described  along  with 
the  equations  of  motion,  a  program  to  calculate  trajectory 
was  possible.  Figure  14  describes  the  process. 

5 .  Validity 

The  program  initially  was  modified  to  print  the  nose 
shapes  every  second  along  the  trajectory  so  that  a  comparison 
between  the  desired  nose  shapes  described  in  Figure  10  and 
the  nose  shapes  described  by  the  equations  derived  in 
Section  IIIA3  could  be  made.  Figure  15  shows  a  family  of 
curves  plotted  from  the  program  output  for  nose  shape  at 
altitudes  closely  corresponding  to  the  altitudes  for  the  family 
of  curves  in  Figure  10.  Although  the  nose  shapes  in  Figure  15 
are  somewhat  more  square  in  appearance  than  those  in  Figure  10, 
the  proper  general  shapes  were  clearly  achieved  with  the  nose 
shape  as  a  function  of  altitude  equations. 

A  goal  of  the  model  was  to  achieve  a  drag  coefficient 
of  about  0.1  and  a  ballistic  coefficit  of  about  2000.  In  the 
simulation,  the  drag  coefficient  ranged  from  0.100  to  0.127 
during  the  flight  with  the  initial  value  of  0.105.  The 
ballistic  coefficient  ranged  from  1392  to  1878  during  the 
flight  with  the  initial  value  of  1798.  The  range  should  be 
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Figure  14.  Flow  Chart  of  Baseline  Trajectory  Program 
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less  than  the  vacuum  trajectory  range  of  101.5  nm;  the  range 
with  drag  was  84.103  nm. 


For  ballistic  entries  at  escape  velocity  and  initial 
angles  of  inclination  greater  than  about  5° ,  the  variation  of 
yelocity  with  altitude  has  been  computed  by  Loh  [Ref.  7]  for 
various  values  of  ballistic  coefficients.  Figure  16  shows  a 
plot  of  the  various  deceleration  curves.  A  similar  curve  can 
be  generated  from  the  simulated  reentry.  The  generated  curve 
should  be  the  same  shape  as  Figure  16  but  since  the  reentry 
velocity  is  only  about  80  percent  of  escape  velocity,  all 
values  of  dimensionless  velocity  will  be  lower. 

In  addition,  the  simulated  reentry  body  undergoes 
ablation  and  hence  changes  the  value  of  the  ballistic 
coefficient.  However,  if  a  nominal  ballistic  coefficient 
of  1500  is  chosen,  the  val".e  of 


WsinO 


1 

SsxnQ 


1500  sin  25‘ 


0.0016 


can  be  used,  where  9  corresponds  to  reentry  angle.  Quantita¬ 
tively,  the  simulated  reentry  curve  should  lie  between  the 
0.001  and  0.01  curves  and  to  the  left  of  the  curves  in 
Figure  16. 

Figure  17  is  a  graph  of  the  Mach  number  versus 
altitude  for  the  model.  Maximum  deceleration  occurs  at 
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25000a 


Figure  17.  Altitude  versus  Mach  Number 


F 


1 


( 
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about  75,000  feet,  and  the  actual  impact  Mach  number  is  4.42. 
Based  on  the  above  parameters  and  the  similarity  in  shape 
between  curves  in  Figures  16  and  17,  the  computer  model  was 
accepted  as  a  valid  representation  of  the  actual  reentry 
trajectory. 


B.  TRAJECTORY  THROUGH  AN  ICE  CLOUD 
1.  Model  Cloud 

In  order  to  evaluate  realistically  the  effects  of 
ice  clouds  on  the  trajectory  of  a  reentry  body,  a  repre¬ 
sentative  ice  cloud  was  needed.  Such  information  as  altitude, 
thickness,  particle  size  and  concentration  was  needed. 

Although  information  such  as  this  is  readily  available  for 
relatively  low  altitude  rain  clouds,  useful  information 
regarding  high  altitude  ice  clouds  is  very  rare.  The  only 
information  found  is  contained  in  Table  5,  reproduced  from 
Heymsfield  [Ref.  8],  in  which  aircraft  measurements  of  five 
different  cirrus  uncinus  clouds  were  made. 
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Table  5 

Summary  of  Cloud  Measurements 
(Reproduced  from  Heymsfield  [Ref.  8],  p.  800) 
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The  table  shows  that  the  clouds  ranged  in  altitude 


! 

from  4.45  to  11  kilometers  (14,600  to  36,100  feet)  and 
varied  in  thickness  from  0.3  to  3.45  kilometers  (980  to 
11,300  feet).  Based  upon  the  above  figures  and  the  advice 
of  Dr.  K.  L.  Davidson,  Meteorology  Department,  Naval 
Postgraduate  School,  the  representative  ice  cloud  was 
selected  to  be  at  an  altitude  of  30,000  feet  with  a  thickness 
of  5,000  feet.  Ice  particle  concentration  at  that  altitude 
nominally  was  chosen  as  0.24  gram  per  cubic  meter  which 
corresponds  to  4.65  x  10  7  slug  per  cubic  foot. 

2 .  Mass  Defect 

The  mass  of  the  reentry  body  ablated  by  collision 
with  the  ice  particles  can  be  expressed  by 


Mass  Defect 


Base  area  x  Cloud  thickness  x  Mass  Loss  Factor 

cos  A 


(58) 


where  the  base  area  of  the  reentry  body  and  the  angle  of  the 
trajectory  with  respect  to  the  horizontal,  A,  are  computed 
at  the  proper  altitude  by  the  computer  program.  The  Mass 
Loss  Factor  represents  the  ratio  of  the  mass  ablated  to  the 
mass  of  the  ice  particles  which  collide  with  the  body.  Some 
information  for  Mass  Loss  Factor  exists  for  multiple  impacts 
without  a  shock  wave  ahead  of  the  reentry  body,  but  none  could 
be  found  which  included  the  shock  wave.  Therefore,  several 
nominal  values  of  Mass  Loss  Factor  were  selected,  and  the 


matrix  of  10,  20,  30,  and  40  was  computed.  A  trajectory 
was  calculated  for  each. 

One  assumption  inherent  in  equation  (58)  is  that  all 
of  the  ice  particles  contained  in  a  tube  representing  the 
bu.  .  area  of  the  reentry  body  passed  at  the  proper  angle 
through  the  cloud  thickness  do  impact  on  the  body.  Such  an 
assumption  was  made  since  this  analysis  was  based  on  a  worst 
case  scenario.  If  only  a  portion  of  the  ice  impacts  on  the 
body,  which  is  the  more  likely  case,  less  ablation  will  occur 
and,  therefore,  less  of  an  effect  will  result. 

3 .  Flow  Chart 

Using  the  cloud  information  contained  in  Table  5 
and  equaticii  (58)  ,  tne  flow  chart  in  Figure  14  was  modified 
to  include  ablation  due  to  ice  erosion.  Figure  18  describes 
the  process. 
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Figure  18. 


Flow  Chart  of  Trajectory  Program  with  Ice  Cloud 
Ablation 
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IV.  RESULTS,  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  RESULTS 

Table  6  compares  the  difference  in  range  of  the  various 
trajectories  modeled  as  functions  of  Mass  Loss  Factor,  which 
is  the  ratio  of  mass  loss  of  ablation  material  due  to  erosion 
to  the  mass  of  ice  particles  impacting.  Also  shown  are  the 
values  of  mass  loss  due  to  ice  erosion  and  final  Mach  number. 
Figure  19  shows  the  same  results  graphically. 


Table  6 

Comparison  of  Trajectories  with  Ice  Erosion 
to  Baseline  Trajectory 


Haaa  Losa  rac tor 

Tina!  Mach  "lumber 

Mas*  Losa  Due  to  lea  (aluaa) 

01 ff trance  In  Range  (feet) 

0 

4.417 

0 

o 

?0 

4.418 

0.0747 

1.54 

20 

4.420 

0.1875 

2.85 

30 

0.2842 

7.99 

40 

4.42'} 

0.3787 

B.  CONCLUSIONS 

The  most  interesting  result  of  the  analysis  is  that  the 
reentry  vehicle  with  ice  erosion  traveled  farther  in  all 
cases  than  did  the  vehicle  without  ice  erosion.  The  reason 
is  that,  despite  the  significant  mass  loss  due  to  erosion 
which  tends  to  decrease  the  ballistic  coefficient,  the  final 


MASS  LOSS  FRCTOR 


nose  shape,  as  seen  in  Figures  10  and  15,  reduces  the  drag 
coefficient  and  tends  to  increase  the  ballistic  coefficient. 

The  result  is  that  the  two  parameters,  mass  loss  and  drag 
coefficient,  are  of  opposite  sign  and  tend  to  cancel. 

Circular  Error  Probability,  CEP,  is  a  constant  multi¬ 
plied  by  the  variance  in  range,  and  for  a  Gaussian  distribution 

a"  =  Z  a.4  (59) 

i 

where  denotes  the  variance  in  range  caused  by  the  various 
factors  that  influence  changes  in  range.  The  assumption  can 
be  made  that  another  ,  that  is,  the  variance  due  to  ice 
erosion,  can  be  added  to  equation  (59)  and  further  that  the 
variance  due  to  ice  erosion  is  of  the  order  of  the  change  in 
range  due  to  ice  erosion.  If  it  also  is  assumed  that  the  CEP 
is  of  the  order  of  the  variance,  then  equation  (59)  can  be 
written 


CEP2  =  Z  CEPi2  =  CEP12  +  A  Range2  (60) 

where  CEP^  is  the  original  CEP  of  the  body  due  to  all 
factors  excluding  ice  erosion  and  A  Range  is  the  difference 
in  range  due  to  ice  erosion.  For  a  C-4  Trident  with  an 
initial  CEP  of  1500  feet  (from  Table  2)  and  a  Mass  Loss 
Factor  of  40  for  a  A  Range  of  15.03  feet,  equation  (60) 
yields  a  new  CEP  of  1500.08  feet. 
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The  assumptions  in  this  thesis  were  made  on  the  basis  of 
the  best  information  available  in  the  literature.  The  analyse 
of  the  problem  based  on  the  assumptions  indicates  that  the 
difference  in  CEP  due  to  ice  erosion  is  very  small  indeed. 
Considering  ice  erosion  as  the  worst  case  for  atmospheric 
water  phenomena,  the  CEP  change  due  to  water  droplet  erosion 
or  due  to  differences  in  humidity  also  can  be  considered  to 
be  very  small. 

C.  RECOMMENDATIONS  FOR  FUTURE  WORK 

1.  The  nose  shape  after  impact  with  ice  particles  is 
obviously  very  important.  Experiments  which  measure  the 
effects  of  multiple  collisions  of  ice  on  nose  shape  are 
required.  Such  experiments  must  include  the  effects  of  a 
bow  shock  around  the  nose  equivalent  to  a  Mach  number  of 
about  10.  Conceptually,  more  ice  should  impact  near  the 
centerline  of  the  vehicle  and  less  toward  the  edges,  tending 
to  make  the  body  mo-e  blunt  and  therefore  tending  to  increase 
the  drag  coefficient. 

2.  This  analysis  should  be  redone  using  data  for  a 
carbon-carbon  nose  and  new  data  for  ice  erosion  found  in 
part  1  above.  The  mass  loss  for  a  carbon-carbon  nose  tip 
would  tend  to  be  far  less  than  for  the  silaceous  nose  tip 
which  is  modeled  in  this  thesis.  Hence,  mass  loss  due  to 
general  erosion  would  tend  to  have  less  effect  than  in  this 
analysis . 


It  may  be,  however,  that  the  mass  loss  due  to  ice  erosion 
for  the  carbon-carbon  nose  is  of  the  same  order  of  magnitude 
as  the  mass  loss  due  to  ice  erosion  for  the  silaceous  nose. 

If  so,  the  effects  of  ice  erosion  may  be  more  important  in 
the  carbon-carbon  nose  than  in  the  silaceous  nose  because 
the  effects  of  mass  loss  due  to  ice  erosion  were  masked  in 
the  silaceous  nose  by  the  substantial  mass  loss  occurring 
throughout  the  reentry. 
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APPENDIX  A 


POLYNOMIAL  REPRESENTATION  OF  NOSE  SHAPES 

The  following  pages  contain  graphs  of  the  various  nose 
shapes.  Each  graph  contains  a  plot  of  the  data  points  taken 
from  the  respective  nose  shape  curves  in  Figure  10  with  the 
aid  of  the  HP-9874A  digitizer  and  a  plot  of  the  polynomials 
which  describe  the  data  points.  The  polynomials  also  are 
displayed  on  the  graph;  the  fifth-order  polynomial  applies 
for  0  <  x  <_  2.4  and  the  first-order  polynomial  applies  for 
x  >  2.4. 

In  the  analysis,  the  polynomial  was  used  to  describe  the 
first  2.4  inches  of  the  nose  shape.  The  linear  approximation 
was  used  to  describe  the  nose  shape  from  2.4  inches  aft  of 
the  nose  tip  to  about  20  inches  aft  of  the  nose  tip.  The 
linear  approximation  to  Figure  A-l  was  used  to  describe  the 
nose  shape  from  about  20  inches  aft  of  the  nose  tip  (where 
the  two  linear  portions  intersect)  to  the  end  of  the  reentry 
body  for  all  nose  shapes. 

The  profiles  for  nose  shapes  with  ice  erosion  use  the 
same  equations  to  describe  the  shape.  The  coefficients  are 
modified  by  substitution  of  a  modified  altitude  into  the 
equations  for  the  coefficients  displayed  in  Appendix  B. 
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Figure  A-3.  Nose  Shape  at  56,250  Feet 


figure  A-4.  Nose  Shape  at  40,100  Feet 
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Figure  A-5.  Nose  Shape  at  25,000  Feet 


APPENDIX  B 


POLYNOMIAL  REPRESENTATION  OF  COEFFICIENTS 

The  following  pages  contain  graphs  of  the  various  nose 
shape  coefficients  defined  in  equations  (49)  through  (52) . 
The  first  six  graphs  correspond  to  the  coefficients  used  in 
equations  (49)  and  (51)  and  apply  for  0  <_  x  <_  2 . 4  inches 
where  x  is  centerline  distance  from  the  nose  tip.  The  last 
two  graphs  correspond  to  coefficients  used  in  equations  (50) 
and  (52)  and  apply  for  x  >  2.4  inches. 
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figure  B-2.  Graph  of  Coefficient  A(4)  versus  Altitude 
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Figure  B-6.  Graph  of  Coefficient  A(0)  versus  Altitude 
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Figure  B-8.  Graph  of  Coefficient  B(0)  versus  Altitude 


APPENDIX  C 


POLYNOMIAL  REPRESENTATION  OF  NOSE  TIP  RECESSION 

The  following  graph  shows  the  evaluation  of  nose  tip 
recession,  the  distance  the  nose  tip  ablates  from  the 
original  location,  as  shown  in  Figure  10.  As  for  the 
coefficients  described  in  Appendix  B,  the  nose  tip  recession 
was  divided  into  two  portions.  A  sixth-order  polynomial 
describes  the  nose  tip  recession  below  90,000  feet  and  a 
linear  model  describes  the  recession  above  90,000  feet. 
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APPENDIX  D 


COMPUTER  PROGRAM  EXAMPLE 

Appendix  D  is  a  sample  of  the  computer  program  used  to 
conduct  the  trajectory  analysis.  The  computations  were 
conducted  on  an  HP-9845A  computer  which  uses  BASIC  computer 
language.  In  general,  the  program  is  self-explanatory  in 
that  the  names  of  variables  are  used  whenever  possible. 

Line  numbers  beginning  with  1XXX  describe  the  initial 
conditions.  The  counter  K  is  zero  until  zero  altitude  is 
reached  and  then  is  equal  to  one.  The  counter  L  is  zero 
during  the  first  iteration  and  allows  the  establishment  of 
initial  volume,  mass,  and  base  area  of  the  reentry  body. 
After  the  first  iteration  of  nose  shape,  L  is  equal  to  1. 

The  counter  Z  is  zero  until  the  altitude  of  the  cloud  is 
reached.  The  value  of  Z  is  changed  to  one  during  the 
iteration  to  find  the  mass  loss  due  to  the  impact  of  the 
ice  in  the  cloud.  After  the  iteration  is  complete,  the 
value  of  Z  is  two.  If  Z  is  set  equal  to  two  initially,  the 
baseline  trajectory  with  no  ice  ablation  will  be  calculated. 

Line  numbers  beginning  with  2XXX  calculate  the  new 
position  of  the  body  as  it  descends.  Line  numbers  beginning 
with  4XXX  calculate  all  of  the  major  parameters  needed  for 
the  calculation  of  pressure  coefficients  and  the  nose  shape 
coefficients.  Line  numbers  beginning  with  5XXX  calculate 
the  actual  nose  shape,  pressure  coefficient,  drag,  volume, 
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etc.  Line  numbers  beginning  with  6XXX  calculate  the  initial 
volume,  mass  and  base  area  and  values  of  drag  coefficient 
and  ballistic  coefficient  during  the  entire  reentry.  Line 
numbers  beginning  with  7XXX  print  the  final  values  obtained 
during  reentry  and  stop  the  program.  Line  numbers  beginning 
with  8 XXX  calculate  the  mass  ablated  during  passage  through 
the  ice  cloud.  Line  numbers  beginning  with  9XXX  calculate 
the  final  range  and  set  the  final  altitude  to  zero  so  that 
calculation  of  other  parameters  at  zero  altitude  is  possible. 
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